The structure, energy, and electronic states of vacancies in Ge nanocrystals 



Kenneth Bayus 

Department of Material Science and Engineering, Cornell University, Ithaca, New York 14850 



o 

as: 
< 

(N 



CO 



0. Paz 

Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019 

S. P. BeckmarQ 

Department of Material Science and Engineering, Iowa State University, Ames, Iowa 50014 

(Dated: August 31, 2010) 

The atomic structure, energy of formation, and electronic states of vacancies in H-passivated Ge 
nanocrystals are studied by density functional theory (DFT) methods. The competition between 
quantum self-purification and the free surface relaxations is investigated. The free surfaces of crystals 
smaller than 2 nm distort the Jahn- Teller relaxation and enhance the reconstruction bonds. This 
increases the energy splitting of the quantum states and reduces the energy of formation to as 
low as 1 eV per defect in the smallest nanocrystals. In crystals larger than 2 nm the observed 
symmetry of the Jahn- Teller distortion matches the symmetry expected for bulk Ge crystals. Near 
the nanocrystal's surface the vacancy is found to have an energy of formation no larger than 0.5 to 
1.4 eV per defect, but a vacancy more than 0.7 nm inside the surface has an energy of formation 
that is the same as in bulk Ge. No evidence of the self-purification effect is observed; the dominant 
effect is the free surface relaxations, which allow for the enhanced reconstruction. From the evidence 
in this paper, it is predicted that for moderate sized Ge nanocrystals a vacancy inside the crystal 
will behave bulk-like and not interact strongly with the surface, except when it is within 0.7 nm of 
the surface. 



PACS numbers: 66.30.Pa, 73.22.-f, 61.46.-w, 61.72.uf 
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I. INTRODUCTION 

Germanium is a particularly attractive material for use 
in semiconducting devices. The charge carriers have a 
high mobility due to their low effective mass and it is 
possible to achieve a high level of n- and p-type dopant 
activation*^ The designers of microelectronics initially 
focused on Si instead of Ge because Ge lacks a native ox- 
ide that can be used as a dielectric. Fortunately this lim- 
itation can be overcome by several techniques that have 
been developed within the last decade: a thin Si over- 
layer can be grown over the Ge so that Si02 can be used 
as the dielectric^ a Ge-oxynitride dielectric layer can be 
grown over the Ge^ or a high-K dielectric crystal, such 
as Zr02, can be used in the device^ These techniques al- 
low for the development of Si/Ge heterostructure devices 
such as mctal-oxide-semiconductor field effect transistors, 
MOSFETs. The heterostructure MOSFET is primarily 
Si so existing fabrication technology can be used, but Ge 
is included as a buried channel between the source and 
drain to allow for high-speed conductivity^— 

Before Ge can become an industrially important mate- 
rial, it must be possible to introduce and control a variety 
of dopant species in the crystal^ Due to the initial chal- 
lenge of finding a suitable dielectric material, there has 
been much less effort in understanding Ge as compared to 
Si. Subsequently much less is known about the control of 
dopants in Ge than in Si. There exists a close relationship 
between impurity diffusion and self-diffusion; therefore, it 
is fundamentally important to understand self-diffusion 



to understand the control of impurity atoms. 

Following Refs. [ll[ and [l2j|, the self-diffusion coeffi- 
cient D (T) is written as a sum of the vacancy (v), in- 
terstitial (i), and direct-exchange (ex) diffusion mecha- 



D (T) = f v (T) C% q (T)D V (T) + 

+ f l (T)Ct q {T)D l {T) + D ex (T) (1) 

where f n (T) are correlation factors, C^ q (T) are the 
equilibrium concentrations of the intrinsic defects, and 
D v (T) are the diffusion coefficients corresponding to 
r\ = v,i,ex. The concentrations and diffusion coefficients 
are expressed in terms of their thermodynamic quanti- 
ties, 
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where AS^ and Ai?7 are the entropy and energy of for- 
mation, AS^ and Ai?^ are the entropy and energy of 
migration, ks is the Boltzmann constant, and K m is a 
constant prefactor. The constant K m is independent of 
temperature and depends on the lattice geometry and the 
vibrational frequencies. The entropy terms include both 
the configurational and the vibrational entropies. The 
energy of formation is determined by the atomic bond- 
ing at the defect site and the energy of migration is de- 
termined by the energy of the saddle point configuration 
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TABLE I: The energies of formation and migration for vacan- 
cies and interstitial defects in Si and Ge. 



Element 


AE) 


AEi 








3.1-3.&S 


3.3i£ 


0.4-1.40±2 


0.45 1 ' 




3.7— 


3.31-3.84i2 


0.43-0.49±2 


0.84±2 


3.492a 
3.5322 


3.27— 








2.3^ 


2.29*2 


0.7^ 


0.5^ 


Ge 


1.7-2.021 
2.42^ 
2.622 
2.56^- 


2.3-4.1^ 
3.55^ 
3.502i 


0.36-0.7-22 





along the minimum energy transition path between sta- 
ble atomic configurations. 

Direct exchange is the slowest diffusion mechanism. 
Straining the lattice to allow the atoms to move past each 
other requires a prohibitively large amount of energy. 
The principal diffusion pathways involve either vacancy 
or interstitial assisted migration. The work presented 
here focuses on these intrinsic point defects. Based on the 
energies listed in TABLE U the self-diffusion coefficient 
in Si is controlled by a self-interstitial kick-out mecha- 
nism at high temperatures (T > 900° C) and vacancy- 
mediated diffusion at lower temperatures.— 

In Ge, vacancy-assisted diffusion is the primary 
mode^ Although the migration barrier for Ge vacan- 
cies and self-interstitials is roughly equivalent, the energy 
of formation for interstitial Ge atoms is approximately 
1 cV greater than vacancies, whereas in Si the energy 
to create vacancies and interstitials is equivalent. In Ge 
self-interstitial atoms will only be formed at very high 
temperatures or after highly energetic processes such as 
irradiation— Therefore, vacancies in Ge are substantially 
more influential than interstitial atoms for assisting dif- 
fusion under thermal equilibrium, as compared to Si. 

The dominance of vacancies-assisted diffusion is ob- 
served experimentally, both for self-diffusion 3 ^ and 
impurity-atom diffusion— The interaction between 
vacancies and impurity atoms is complicated. It is be- 
lieved that vacancies and impurities form mobile defect 
pairsi 28 i 35 ' 36 This defect pair can become pinned when a 
second impurity, such as C, joins the complex— In addi- 
tion to forming complexes, the vacancies and impurities 
often carry a charge i 23 i 28 ' 35 It is likely that an isolated va- 
cancy in bulk Ge is charged —2— In the work presented 
here only isolated, charge neutral, impurities are investi- 
gated, which is consistent with the nanoscale context of 
this study. 

To improve the engineering control of material prop- 
erties and increase device efficiency, it is desirable to 
move from bulk to nanoscale structures. There are many 
examples of situations where nanostructured Ge offers 
benefits. The use of Ge nanocrystals as the floating 
gate of MOS memory devices results in a dramatic shift 
in the threshold-voltage, improved switching character- 



istics, and decreased leakage current— ^ Ge films with 
nanostructured surfaces offer the ability to tune the opti- 
cal properties of thin films— Ge nanowires are considered 
for use as MOSFETs4Mi 

Ultimately the nanostructures used to create devices 
need to be tailored by controlling their size, surfaces, 
and dopants. With respect to introducing dopants, the 
electronic properties of nanostructures are believed to 
be sensitive to the relative position of the impurities in 
the structure. The mean free path of charge carriers 
within nanowires depends strongly on the radial dopant 
profile— This will influence the conductivity. In addition 
to the challenge of selectively incorporating the dopant 
atoms into the nanostructure, the impurity distribution 
must be maintained for the lifetime of the device. 

At the quantum scale the primary difference between 
a bulk crystal and a nanocrystal is the interaction of 
the wave function with the surfaces. As the size of a 
structure decreases, the crystal's translational symmetry 
ceases to be meaningful. The electronic band-structure 
that is nominally a function of the quantum number k is 
projected onto the T point in the center of the Brillouin 
zone. The crystal's energy bands become discrete quan- 
tum energy states. Whereas in bulk the wave function 
is distributed across the entire crystal as Bloch waves, 
Mk (r) e 4k ' r , in nanocrystals the wave function is confined 
by the surfaces. The size of the nanostructure directly 
impacts the energy states, analogous to the elementary 
particle-in-a-box problem. Consider for example a [110] 
Ge nanowirc. When the wire diameter is sufficiently 
small the crystal's translational symmetry is only mean- 
ingful in the [110] direction and the bulk Ge states are 
projected along the k = [110] direction in k-space. This 
projection transforms Ge from an indirect to direct band 
gap material— The confinement is predicted to distort 
the shape of the energy dispersion for wires with diam- 
eters as large as 2 nm. The energy bands of nanowires 
with diameter greater than 2 nm are found to undergo a 
rigid shift, even for wires as large as 5 nm— 

Fundamentally, there are two effects that differenti- 
ate the behavior of defects in nanostructures from bulk: 
quantum confinement and free surfaces. Dalpian^i claims 
that the confinement of the defect's wave function results 
in the so-called self-purification effect that increases the 
defect's A£^. In the case of dopant species, this increase 
hinders the incorporation of dopant atoms into the nanos- 
tructures. This is a controversial subject and worthwhile 
investigating—"— In the present calculations evidence of 
self-purification will be sought. 

The free surfaces allow the nanostructure to expand or 
contract to reduce the strain energy surrounding the de- 
fect. From an energetics perspective the self-purification 
effect and the free surfaces compete with one another. 
The self-purification increases the energy and the free 
surfaces decrease the energy. From a kinetics perspective 
they complement one another because it is likely that 
the surfaces will getter impurities out of the nanostruc- 
ture. In the case of Si nanocrystals it is observed that 
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the relative energy to introduce vacancies decreases as 
the nanocrystal's size decreases. This indicates that en- 
ergetically the free surfaces dominate the self-purification 
effect i& As the vacancy is moved toward the surface the 
energy further decreases and when the vacancy is within 
0.6 A of the surface it becomes unstable and is sponta- 
neously moved to the surface of the crysta h 48 ' 49 

In this paper the structure and energies of vacancies 
in Ge nanocrystals are examined as a function of the 
nanocrystal's size and the position of the vacancy in the 
crystal. Because the energies, A£?y, depend on the size 
of the crystal and the position within the crystal, the 
concentrations and diffusivities, from EQ. [5] and [31 also 
depend on the size and position. Self-diffusion within 
nanostructures is not a simple matter that can be easily 
described by a single coefficient. The work here is a first 
step toward building a comprehensive model. 

Following this introduction, in Sec. HU the methods 
used will be presented, including a discussion of the com- 
putational approach and the details of the nanocrystals' 
morphology. The results from these calculations will be 
presented and discussed in Sec. IIIII A concluding sum- 
mary will be presented in Sec. IIVI 



II. METHODS 

A. Computational approach 

The calculations are performed within the frame- 
work of the density- functional theory^ (DFT) , using the 
local-density approximation^ 1 - (LDA) for the exchange- 
correlation functional, as it is implemented in the Siesta 
computational packaged The electrons in the core 
atomic region are substituted by norm-conserving pseu- 
dopotcntials of the Troullier-Martins type,— and the va- 
lence charge is represented by a set of atom-centered ba- 
sis functions. In Siesta these functions correspond to 
numerical atomic orbitals of strictly finite range, a par- 
ticular choice that is specially suited to treat isolated 
systems. 

All calculations are carried out using a double-^ plus 
polarization orbitals (DZP) basis set. The cutoff radii of 
the basis functions are optimized for bulk germanium in 
the diamond structure, following the method proposed 
in Ref. [54|, where a fictitious external pressure of 0.2 
GPa is employed on the free atom. This basis size and 
radius lengths are proven to give a good balance between 
the computational accuracy and cost. A theoretical lat- 
tice parameter of do = 5.64 A and a bulk modulus of 
Bo = 80.0 GPa, are obtained from the fitting to the 
Murnaghan equation of state.— Both values are in good 
agreement with the structural and elastic properties from 
experiments^ (<io = 5.66 A and Bq = 75.8 GPa), given 
the fact that the LDA tends to underestimate lattice con- 
stants by a 1-3%, but also to overestimate bulk moduli 
with errors ranging from 5 up to 20% J% 

The theoretical method employs periodic boundary 



TABLE II: The local bond lengths at a vacancy site in crystals 
with diameters 1.02 and 2.20 nm. The segment labels refer- 
ence the tetrahedral structure in FIG. [2] Using the DFT-LDA 
the theoretical bond length in bulk Ge is 2.44 A, which corre- 
sponds to segment lengths of 3.99 A before atomic relaxation. 



Segment D = 1.02 nm D = 2.20 nm 

AB 2.55 2.90 

AC 3.57 3.47 

AD 3.57 3.47 

BC 3.79 3.46 

BD 3.79 3.46 

CD 2.50 2.90 



conditions. The nanocrystals are placed inside the su- 
percell surrounded by a buffer of empty space. The size 
of the atomic clusters ranges from 44 to 244 Ge atoms 
and the vacuum region is chosen to be large enough as 
to avoid any interaction between their periodic replicas. 
A kinetic energy cutoff of 250 Rydbcrg is chosen for 
the real-space integrations involving the Hartrce and the 
exchange-correlation contributions to the self-consistent 
potential. In this respect, a stringent criterion is em- 
ployed in the convergence of the density matrix and total 
energy. All atomic coordinates are then relaxed accord- 
ing to a conjugate-gradient minimization algorithm, until 
the maximum residual forces are below 0.02 eV/A. 



B. Nanocrystal morphology 

The nanocrystal geometries used in these studies 
are hydrogen passivated, bond-centered crystals. Ex- 
perimental nanostructures frequently have amorphous, 
glassy, or polymeric coatings that result from the method 
of crystal growth. It is possible to treat the surfaces to 
reduce or remove these, although it is uncommon exper- 
imentally to work with bare surfaces. The nanocrystals 
investigated here have their surfaces passivated with a 
extremely "soft" H pseudopotential. Surface passivation 
removes surface states and allows the competition be- 
tween the self-purification effect and the free surfaces to 
be studied without considering the complicated surface 
chemistry. 

Surface Ge atoms are identified and the dangling bonds 
of the Ge atoms are capped with H. Any Ge atom that is 
found to have three dangling bonds is replaced by a single 
H atom. The resulting nanocrystals are shown in FIG. 
[TJ The surface morphologies are examined and crystals 
that are highly faceted are excluded. Only nanocrystals 
with near spherical geometry are studied. 



FIG. 1: (color online) Ge nanocrystal morphology for (a) Ge44ri42, (b) Gei3oHg8, and (c) Ge244Hi58. The corresponding sizes 
are 1.02, 1.70 and 2.20 nm respectively. A central vacancy is depicted as a hatched atom surrounded by four missing bonds 
(broken lines), all in red. Ge atoms are colored in blue; saturating H atoms are in light grey. 




FIG. 2: The reference geometry of the atomic structure at 
the vacancy site. The dashed circle is the vacancy and the 
solid circles are the nearest neighbor Ge atoms that form a 
tetrahedron. The tetrahedron, without atomic motion has 
Td symmetry. The calculated bond lengths for various sized 
nanocrystals are shown in TABLE ITT1 



III. RESULTS AND DISCUSSION 
A. Atomic structure and defect states 



An undistorted vacancy has Td symmetry. There are 
three degenerate states in the gap, belonging to the t-z 
representation, associated with this structure. These are 
shown in the middle column of FIG. |3] There are two 
electrons localized at the defect so the states are par- 
tially occupied. It is the partial occupancy of the degen- 
erate energy levels that allows the defect to undergo a 
spontaneous symmetry breaking that reduces the degen- 
eracy and lowers the electronic energy. In bulk crystals 
the Jahn- Teller distortion produces a D 2 d symmetrized 
structure with the fully-occupied state belonging to the 
&2 representation and the doubly-degenerate, empty state 
having the e representation i 48 i 58 i 59 

Here, as in the case of a vacancy in a Si nanocrystal^ 8 - 
the symmetry of the structure approaches D 2 d but due to 
the surfaces there is additional distortion. In the case of 
the 1.02 nm crystal, the symmetry of the vacancy struc- 
ture is C s . For the 2.20 nm crystal the symmetry is 
essentially D 2 d but a minuscule 0.01 nm distortion of the 
bonds lowers the symmetry, i.e., if AC=AD=BC=BD 
then the symmetry would be D 2 d- From TABLE [XT] it 
is determined that the Jahn- Teller distortion in the 1.02 
nm crystal is approximately 10% larger than that in the 
2.20 nm crystal and as a result the defect states undergo 
a larger energy split, over 2.7 eV, which almost pushes 
the states out of the gap. In the 2.20 nm crystal the 
splitting is much smaller, around 0.75 eV. 



The local atomic structure at the vacancy site is di- 
rectly related to the electronic states introduced to the 
gap from the broken bonds. The atomic structure of the 
vacancy in the 1.02 and 2.20 nm nanocrystals is given in 
TABLEAU using as reference the ABCD indices shown in 
FIG. O The associated electronic states are diagrammed 
in FIG. El The left column of FIG. [3] shows the states 
for a nanocrystal with no vacancy. The band gap for the 
1.02 nm crystal is 3.13 eV and the gap for the 2.20 nm 
crystal is 2.0 eV. 



B. Crystal size 

The energy of the fully optimized vacancy structures in 
the different sized nanocrystals is calculated. Subtract- 
ing this energy from the energy of the perfect nanocrys- 
tals yields the energy of formation for a vacancy plus 
the chemical potential for Ge, i.e., the energy to remove 
a Ge atom from the system. The chemical potential is 
variable and depends on the local chemical environment. 



Vacancy Vacancy 
Perfect Crystal No Relaxation Jahn-Teller 
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FIG. 3: The top frame shows the energy levels for a crystal 
with 1.02 nm diameter and the bottom frame shows a 2.20 nm 
crystal. The red dashed line is the Fermi energy. The band 
alignment between the diagrams is arbitrary. The crystals 
in the left column are perfect and contain no vacancy. The 
crystals in the center column have vacancies, but the struc- 
ture has not been allowed to relax. The states at the Fermi 
level are three-fold degenerate. The right column shows the 
states after the Jahn-Teller distortion is completed, which al- 
lows the atomic structure to break the symmetry and lift the 
degeneracy of the gap state. 



By assuming that all the nanocrystals are located in the 
same environment and have the same chemical potential 
it is possible to compare the relative energy of formation 
for vacancies in different sized nanocrystals. It is known 
that as a nanocrystal's diameter approaches infinity the 
energy of formation approaches that found in bulk Gc. 
Using this energy limit the calculated energy of forma- 
tion versus crystal diameter is plotted in FIG.0|a). It is 




Diameter (nm) 



FIG. 4: The relative energy of formation for a vacancy in 
different sized nanocrystals is plotted in frame (a). In frame 
(b) band gap is plotted versus the crystal size. The solid lines 
show EQ.[3]with the coefficients fitted to the data. 



assumed that the size dependence goes as 

E{D) = ^ +1 (4) 

where a, /3, and 7 are fitting coefficients. In the limit 
that the diameter, D, goes to infinity, the energy equals 
7. The coefficients are determined to be a = —1.1395 and 
(3 = 6.2574. The energy zero is shifted so that 7 = 2.0 
eV, which is taken from the energies reported in TABLE 
H The quality of this fit appears good. It is observed that 
the energy of formation is near the bulk value for crystals 
as small as 2.0 nm. This is surprising because quantum 
confinement continues to strongly influence the band gap 
for crystals with a similar size, as shown in FIG. |4|b). 



C. Distance from crystal center 

To determine the influence of a vacancy's position on 
its energy a 2.20 nm crystal (GC244H158) is examined with 
a vacancy at various locations within it. The calculated 
energies are shown in FIG. [5j The configuration where 
the vacancy is adjacent to the center of the crystal [FIG. 
[He)] is defined as the zero. Near the center of the crystal 
there is little change in the energy, but once the vacancy 
is within 0.7 nm of the surface it begins to drop substan- 
tially The last stable vacancy site is 0.3 nm from the 
surface. The energy of a vacancy at this site is a full 1.2 
eV less than a vacancy near the center. From the results 
in Sec. MI Bl it is known that the energy of formation in 
the center of the 2.20 nm crystal is almost that observed 
in bulk or slightly smaller. Using the energies in TABLE 
U it is deduced that the energy of formation for the va- 
cancies near the surface can be no larger than 0.5 to 1.4 
eV. 
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Vacancy energy versus position 

22 A Ge nanocrystal 




2 4 6 

Distance from center (A) 

FIG. 5: The relative energy of a vacancy in a 2.20 nm Ge 
nanocrystal (Ge244Hi5g). The zero of the energy scale is de- 
fined to be the inner most atomic site. 



IV. CONCLUSION 

A vacancy in a Ge nanocrystal undergoes a Jahn- Teller 
distortion. The Td symmetrized broken bonds located 
at the vacancy introduce a set of three-fold degenerate, 
partially-occupied, states in the gap. When these dan- 
gling bonds reconstruct the defect symmetry is lowered. 
This reduces the degeneracy of the defect states by split- 
ting them into a lower-energy, fully-occupied state and 
two higher-energy, degenerate empty states. In a bulk 
crystal it is known that the symmetry of the vacancy site 
is -D2d,— but in the nanocrystals the surfaces introduce 
additional distortion. For the smallest crystal the surface 
influence is great; the defect symmetry is C s and the en- 
ergy splitting is approximately 2.7 cV. This results in a 
dramatic reduction in energy of formation. In the 2.20 
nm crystal the defect almost has the Did symmetry that 
is found in bulk. The energy splitting is also smaller than 
that in the 1.02 nm crystal, only around 0.75 eV; there- 
fore, the energy reduction due to the bond reconstruction 
is lower and the energy of formation is larger in the 2.20 
nm crystal. This is consistent with the calculated predic- 
tion that the energy of formation will approach the bulk 
value for nanocrystals larger than 2.0 nm. 

The band gap of the crystal continues to change greatly 



even when the diameter is as large as 2 nm. It is deduced 
that although quantum confinement continues to impact 
the energy levels in the crystal, the primary influence on 
the vacancy is the surface's ability to enhance the internal 
structural relaxation. It is concluded that in this example 
the quantum self-purification effect plays a small role if 
any. A similar observation has been made for vacancies 
in Si4£ It is hypothesized that this is due to the defect's 
wave function being highly localized at the reconstruction 
bonds. 

Finally, it is determined that vacancies placed within 
0.7 nm of the surface are spontaneously removed. Sur- 
prisingly vacancies in the interior of the crystal are stable 
and do not appear to be drawn toward the exterior. An 
additional consequence is that if a surface were to act as 
a vacancy source, the vacancies produced from the sur- 
face are unlikely to penetrate deeply into the nanocrys- 
tal. The system studied here has H-passivated surfaces, 
which allows for large relaxations. Experimental crystals 
that have surface reconstructions or polymer coatings will 
have more rigidity and the influence of the surfaces will 
be further muted inside the crystal. 

The picture that emerges from this work is that mod- 
erate sized crystals will have an interior where vacancies 
behave bulk-like and a thin exterior surface region where 
the surface effects will dominate. Assuming that the 
properties of the self-interstitial defect are not strongly 
modified by the surfaces, then the evidence in this pa- 
per predicts that the self-diffusion in the interior of Ge 
nanocrystals will not be substantially different from that 
observed in bulk. However, recent experiments indicate 
that Ge surfaces are not sinks for interstitial atoms, but 
instead reflect the interstitial Ge back into the crystal^ 
If this observation holds within the nano-regime then 
it is possible that the large surface to volume ratio in 
nanostructures will magnify the impact of the interstitial- 
assisted diffusion. 
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